Study of J/<0 



M. Ablikim 1 , J. Z. Bai 1 , Y. Ban 10 , J. G. Bian 1 , D. V. Bugg 19 , X. Cai\ 
J. F. Chang 1 , H. F. Chen 16 , H. S. Chen 1 , H. X. Chen 1 , J. C. Chen 1 , Jin Chen 1 , 
Jun Chen 6 , M. L. Chen 1 , Y. B. Chen 1 , S. P. Chi 2 , Y. P. Chu 1 , X. Z. Cui 1 , 
H. L. Dai 1 , Y. S. Dai 18 , Z. Y. Deng 1 , L. Y. Dong 1 , S. X. Du 1 , Z. Z. Du 1 , J. Fang 1 , 
S. S. Fang 2 , C. D. Fu 1 , H. Y. Fu 1 , C. S. Gao 1 , Y. N. Gao 14 , M. Y. Gong 1 , 
W. X. Gong 1 , S. D. Gu 1 , Y. N. Guo 1 , Y. Q. Guo 1 , Z. J. Guo 15 , F. A. Harris 15 , 
K. L. He 1 , M. He 11 , X. He 1 , Y. K. Heng 1 , H. M. Hu 1 , T. Hu 1 , G. S. Huang 11 ' , 
L. Huang 6 , X. P. Huang 1 , X. B. Ji 1 , Q. Y. Jia 10 , C. H. Jiang 1 , X. S. Jiang 1 , 
D. P. Jin 1 , S. Jin 1 , Y. Jin 1 , Y. F. Lai 1 , F. Li 1 , G. Li 1 , H. H. Li 1 , J. Li 1 , J. C. Li 1 , 

Q. J. Li 1 , R. B. Li 1 , R. Y. Li 1 , S. M. Li 1 , W. G. Li 1 , X. L. Li 7 , X. Q. Li 9 , 
X. S. Li 14 , Y. F. Liang 13 , H. B. Liao 5 , C. X. Liu 1 , F. Liu 5 , Fang Liu 16 , H. M. Liu 1 , 
J. B. Liu 1 , J. P. Liu 17 , R. G. Liu 1 , Z. A. Liu 1 , Z. X. Liu 1 , F. Lu 1 , G. R. Lu 4 , 
_i ' J. G. Lu 1 , C. L. Luo 8 , X. L. Luo 1 , F. C. Ma 7 , J. M. Ma 1 , L. L. Ma 11 , Q. M. Ma 1 , 

X. Y. Ma 1 , Z. P. Mao 1 , X. H. Mo 1 , J. Nie 1 , Z. D. Nie 1 , S. L. Olsen 15 , 
H. P. Peng 16 , N. D. Qi 1 , C. D. Qian 12 , H. Qin 8 , T. N. Ruan 16 J. F. Qiu 1 , 
Z. Y. Ren 1 , G. Rong 1 , L. Y. Shan 1 , L. Shang 1 , D. L. Shen 1 , X. Y. Shen 1 , 
H. Y. Sheng 1 , F. Shi 1 , X. Shi 10 , H. S. Sun 1 , S. S. Sun 16 , Y. Z. Sun 1 , Z. J. Sun 1 , 
X. Tang 1 , N. Tao 16 , Y. R. Tian 14 , G. L. Tong 1 , G. S. Varner 15 , D. Y. Wang 1 , 
O ; J- Z. Wang 1 , K. Wang 16 , L. Wang 1 , L. S. Wang 1 , M. Wang 1 , P. Wang 1 , 

P. L. Wang 1 , S. Z. Wang 1 , W. F. Wang 1 , Y. F. Wang 1 , Zhe Wang 1 , Z. Wang 1 , 
Zheng Wang 1 , Z. Y. Wang 1 , C. L. Wei 1 , D. H. Wei 3 , N. Wu 1 , Y. M. Wu 1 , 
J^; X. M. Xia 1 , X. X. Xie 1 , B. Xin 7 , G. F. Xu 1 , H. Xu 1 , Y. Xu 1 , S. T. Xue 1 , 

M. L. Yan 16 , F. Yang 9 , H. X. Yang 1 , J. Yang 16 , S. D. Yang 1 , Y. X. Yang 3 , M. Ye 1 , 
M. H. Ye 2 , Y. X. Ye 16 , L. H. Yi 6 , Z. Y. Yi 1 , C. S. Yu 1 , G. W. Yu 1 , C. Z. Yuan 1 , 
J. M. Yuan 1 , Y. Yuan 1 , Q. Yue 1 , S. L. Zang 1 , Yu. Zeng 1 ^. Zeng 6 , B. X. Zhang 1 , 
B. Y. Zhang 1 , C. C. Zhang 1 , D. H. Zhang 1 , H. Y. Zhang 1 , J. Zhang 1 , 
J. Y. Zhang 1 , J. W. Zhang 1 , L. S. Zhang 1 , Q. J. Zhang 1 , S. Q. Zhang 1 , 
X. M. Zhang 1 , X. Y. Zhang 11 , Y. J. Zhang 10 , Y. Y. Zhang 1 , Yiyun Zhang 13 , 
Z. P. Zhang 16 , Z. Q. Zhang 4 , D. X. Zhao 1 , J. B. Zhao 1 , J. W. Zhao 1 , M. G. Zhao 9 , 
P. P. Zhao 1 , W. R. Zhao 1 , X. J. Zhao 1 , Y. B. Zhao 1 , Z. G. Zhao 1 *, H. Q. Zheng 10 , 
J. P. Zheng 1 , L. S. Zheng 1 , Z. P. Zheng 1 , X. C. Zhong 1 , B. Q. Zhou 1 , 
G. M. Zhou 1 , L. Zhou 1 , N. F. Zhou 1 , K. J. Zhu 1 , Q. M. Zhu 1 , Y. C. Zhu 1 , 
Y. S. Zhu 1 , Yingchun Zhu 1 , Z. A. Zhu 1 , B. A. Zhuang 1 , B. S. Zou 1 . 

(BES Collaboration) 

1 Institute of High Energy Physics, Beijing 100039, People's Republic of China 

2 China Center for Advanced Science and Technology (CC AST), Beijing 100080, 

People's Republic of China 

3 Guangxi Normal University, Guilin 541004, People's Republic of China 

4 Henan Normal University, Xinxiang 453002, People's Republic of China 
5 Huazhong Normal University, Wuhan 430079, People's Republic of China 

6 Hunan University, Changsha 410082, People's Republic of China 
7 Liaoning University, Shenyang 110036, People's Republic of China 



J3 



Preprint submitted to Elsevier Preprint 



7 February 2008 



Nanjing Normal University, Nanjing 210097, People's Republic of China 

9 Nankai University, Tianjin 300071, People's Republic of China 

10 Peking University, Beijing 100871, People's Republic of China 
11 Shandong University, Jinan 250100, People's Republic of China 

12 Shanghai Jiaotong University, Shanghai 200030, People's Republic of China 

13 Sichuan University, Chengdu 610064, People's Republic of China 

14 Tsinghua University, Beijing 100084, People's Republic of China 

15 University of Hawaii, Honolulu, Hawaii 96822, USA 
16 University of Science and Technology of China, Hefei 230026, People's Republic 

of China 

17 Wuhan University, Wuhan 430072, People's Republic of China 
18 Zhejiang University, Hangzhou 310028, People's Republic of China 
19 Queen Mary, University of London, London El 4NS, UK 

* Current address: University of Michigan, Ann Arbor, Michigan, 48109, USA 
t Current address: Purdue University, West Lafayette, Indiana 47907, USA. 



Abstract 

New data are presented on J/ip — ► ujK + K~ from a sample of 58M J/ijj events in the 
upgraded BES II detector at the BEPC. There is a conspicuous signal for /o(1710) — ► 
K + K~ and a peak at higher mass which may be fitted with /2(2150) — > KK. 
From a combined analysis with ojti + it~ data, the branching ratio BR(fo(1710) — > 
7T7r)/Bi2(/o(1710) -» KK) is < 0.11 at the 95% confidence level. 



PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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In a recent publication, we have presented new data on J/ip — > ujtt + 7t~ [1] 
from a sample of 58M J/^ events taken in the Beijing Spectrometer (BES) 
detector at the Beijing Electron Positron Collider. Here we report data on 
J ftp — * uoK + K~ . Earlier data on this channel with lower statistics have been 
published by Mark I [2], DM2 [3] and Mark III [4]. 

The BES II detector is a large solid-angle magnetic spectrometer that is de- 
scribed in detail in Ref . [5] . Charged particles are measured in a vertex cham- 
ber and Main Drift Chamber (MDC); these are surrounded by a solenoidal 
magnet providing a nearly uniform field of 0.4T. Photons are detected in a 
Barrel Shower Counter (BSC) made of gas proportional tubes interleaved with 
12 radiation lengths of lead sheets. A time-of-flight (TOF) hodoscope immedi- 
ately outside the MDC provides separation between pions, kaons, and protons. 
The time resolution of the TOF measurement is 180 ps. Further separation is 
obtained using dE/dx in the MDC. 

The point of closest approach of a charged track to the beam is required to be 
within 2 cm of the beam axis and within 20 cm of the centre of the interaction 
region along the beam axis. Both photons are required to be isolated from 
charged tracks by demanding an angle > 8° to the nearest charged track. Any 
photon with an energy deposit < 30 MeV in the shower counter is rejected. 
All particles are required to lie well within the acceptance of the detector, 
with charged tracks having laboratory polar angles 6 satisfying | cos#| < 0.84 
and with transverse momenta > 60 MeV/c. 

The uj is observed decaying to 7r + 7r~7r°, so events are selected initially by 
demanding two photons and four charged tracks with total charge zero. If 
there are more than two photons, all are tried; an extra photon can arise from 
interactions of charged particles in the detector. Kaons can be identified up to 
momenta of 800 MeV/c by TOF and dE/dx measurements. The two slowest 
particles always have energies < 800 MeV. The first step is to identify one 
kaon and one pion using TOF and dE/dx. The other two tracks often have 
momenta too high to be identified by TOF and dE/dx, so a four-constraint 
kinematic fit is made for the K + fT~7r + 7r~77 hypothesis. The kinematic fit 
requires x 2 (^ + ^~' K+ ' K ~n) < 40. 

The 7T° is selected by requiring |M 77 — M n o\ < 0.020 GeV/c 2 ; the 7r° mass 
resolution is ~ 15 MeV/c 2 . The resulting 7T + 7r 7T° mass distribution is shown 
in Fig. 1. The uj signal is then selected requiring \M^+^-^o—M w \ < 40 MeV/c 2 . 
The background is fitted by a second order polynomial in M(7r + 7r~7r°). A 
background of (22.9 ± 2.0)% is estimated from uj sidebands, defined by 80 < 
\M n + n - n o — M u \ < 160 MeV/c 2 ; the error allows for small variations when the 
location and width of the sidebins are changed. 

For a given uj momentum, the mass of the accompanying KK pair is unique. 
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The decay angles of irir and KK in the lab frame are very different except 
near or 180°. There, the backward 7r or K differ strongly in momentum and 
are easily distinguished by momentum, TOF, and dE/dx. As a result, there 
is a clean separation between co>7r + 7r~ and loK + K~ . 

Most background originates from K + K~ -k + -k~ it® . The other sources of back- 
ground are Kg in final states K s K ± 7i T 7i° and KgK ir^'y. Most K$ events are 
rejected as follows. ll^{KgK ± 'K T ^) < x 2 (K + K~ 7 r + 7i~<yy) or x 2 {K s K ± -k t ^) < 
X 2 {K + K~7T + 7T~'-y'y), events are discarded if any Ktttttt combination has M(7r + 7r~) 
in the interval 497 ±25 MeV/c 2 and r xy > 3 mm; here r xy is the distance from 
the beam axis to the 7r + 7r~ vertex. This avoids rejecting too many signal 
events; surviving K$ background is too small to be visible. The beam spot 
has a o x of 0.6 mm, and the resolution of the second vertex is 1.2 mm in 
xy. After the background subtraction, there are 3438 signal events. From the 
Monte Carlo simulation, the average detection efficiency is 4.0%. 




Fig. 1. The u peak from the selection described in the text; background is estimated 
from the lower curve. 

Fig. 2(a) shows the experimental Dalitz plot, and Figs. 2(c) and (d) show pro- 
jections on to masses of K + K~ and lvK; the shaded area indicates background 
events from the sideband estimation. 

The channels fitted to the data are: 



^cu/o(980) 

-w/ (1710) 

^/ 2 (1270) 
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Fig. 2. (a) and (b): Measured and fitted Dalitz plots for ujK + K~; (c) and (d) are 
projections on to K + K~ and uoK mass. In the latter, histograms show the maximum 
likelihood fit; the shaded region indicates the background estimated from sidebins; 
the dashed curve in (d) shows the magnitude of the ifi(1400) contribution and a 
Kuj contribution at 1945 MeV/c 2 ; (e) and (f) show mass projections of /o and ji 
contributions to K + K~ . The dashed curve of (e) shows the a — > K + K~ S-wave 
contribution. 
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^/2(1525) or w/ 2 (1565) 

->a;/ 2 (2150) 

-> #i ( 1400) X 

->/ri(i950)iir. 

Amplitudes are fitted to relativistic tensor expressions documented in Ref. [6]. 
For spin in KK, two transitions from J/ip are allowed with orbital angular 
momenta £ = and 2 in the production process. For spin 2, there are five 
amplitudes: one with £ = 0, three with £ = 2 and one with £ = 4. In fitting 
these, Blatt-Weisskopf centrifugal barrier factors are included with a radius of 
0.8 fm, though results are insensitive to this choice. In the amplitude analysis, 
information from the u — > 7r + 7r~7r° decay is included in the tensor expressions. 

The polarisation vector of the uj lies along the normal to its decay plane. The 
correlation between this polarisation vector, the production plane, and the 
decay of the fj to K + K~ is sensitive to the spin of fj and also to the helicity 
amplitudes for its production. This correlation cannot readily be displayed, 
since it depends on five angles; however, tests with different J p demonstrate 
the sensitivity to quantum numbers. 

Fig. 2(b) shows the Dalitz plot from the log likelihood fit described below. 
Histograms on Figs. 2(e) and (f) show projections of / and f 2 contributions 
to this fit. 

The ujit + it~ data of Ref. [1] determine all helicity amplitudes for production 
of / 2 (1270) well. In fitting present data, the relative magnitudes of these am- 
plitudes are fixed to values from ujtctt. Contributions from / (980) are likewise 
fixed from the signal observed in cj7r + 7r~; its branching ratio KK/tttt is taken 
from the Flatte formula fitted to J/x/j — > 07r + 7r~ and <j)K + K~ [7], where there 
are conspicuous / (980) signals. Phases for / 2 (1270) and / (980) amplitudes 
are fitted freely, since they arise from multiple scattering, which is different in 
KK and tctt final states. 

For other components, there is a general problem in isolating / from possible 
f 2 for two reasons. Firstly, five 2 + amplitudes can simulate two + amplitudes 
closely; amplitudes with J p = 2 + may be identified if they give rise to decay 
angular distributions which are non-isotropic. Secondly, fitted 2 + amplitudes 
can fluctuate for angles outside the acceptance. For high K + K~ mass above 2 
GeV/c 2 , this latter problem is somewhat reduced, because the £ = 4 amplitude 
is suppressed by the strong centrifugal barrier for production. 

We use a to denote a broad K + K~ S-wave contribution. We find that it peaks 
towards the lower KK masses as shown by the dashed curve of Fig. 2(e). 
However, the dependence on mass above 1 GeV is somewhat uncertain. Many 
alternative fits have been tried with similar results. A component peaking 
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towards threshold is required; without it, the fit to the KK mass distribution 
of Fig. 2(c) is bad. We have therefore tried parametrisations using the a pole 
of Ref. [1], and a coupling constant of the form G\ + G 2 s or G\ + G 2 /s. 
The optimum fit requires a slightly more rapid fall with s than the a pole, 
in order to fit four points at the lowest KK masses. However, we regard 
this as unphysical and therefore eventually choose to use the a pole of Ref. 
[1] unchanged, with G 2 = 0. Note that there is a substantial constructive 
interference in present data between / (980) and a amplitudes at masses close 
to threshold. 

A dominant feature is /o(1710); the present data are consistent with earlier 
studies which identify J = [8,9]. They are also consistent with the ab- 
sence of any significant J = 2 contribution. The fitted / (1710) optimises at 
M = 1738 ± 30 MeV/c 2 , T = 125 ± 20 MeV/c 2 . The error in the mass is 
mostly systematic, and arises from uncertainty in the a amplitude with which 
/o(1710) interferes; the error in T is mostly statistical, but includes allowance 
for interference with the remaining + amplitude. Earlier BES II data on 
Jji> -> 1 K + K~ and ^K° S K° S gave M = 1740 ± 4(stat)i^(syst) MeV/c 2 and 
r = 166±| t\l MeV/c 2 [8]. 

A fit to the 1738 MeV/c 2 peak with spin 2 uses five amplitudes and gives log 
likelihood worse than spin by only 15; the fit is shown in Fig. 3. However, 
the fit with spin uses only two production amplitudes with £ = and 2. The 
fit with spin requires an £ = amplitude which is completely dominant over 
£ = 2. However, for spin 2 the £ = 2 amplitudes dominate over £ = 0. The 
phase space available in the process J/ip — > cj/j(1710) is rather limited, and 
the £ = 2 and 4 centrifugal barriers for the production process should suppress 
those amplitudes strongly. If the £ = 2 and 4 amplitudes are removed, spin 
gives a fit better than spin 2 by 90 in log likelihood. 

This pattern of behaviour is symptomatic of what is required for spin 2 to 
simulate spin 0. The spin 2 amplitude with £ = has a unique dependence 
on angles; it contains a distinctive term 3cos 2 ai<- — 1, where ax is the decay 
angle of the K + in the resonance rest frame, with respect to the direction of 
the recoil oo. Simulation of spin requires large J = 2 £ = 2 and 4 amplitudes 
to produce compensating terms in sin 2 ax- Although this is suspicious, the 
J = 2 possibility cannot be ruled out from present data. 

We discuss next the branching ratio of / (1710) between KK and rnr, us- 
ing information from J/ip — > wk + ty^ [1], where statistics of ~ 40i^ events 
are available. In those data, there is no definite evidence for the presence of 
/o(1710) ; if its mass is scanned, there is no optimum around 1710 MeV/c 2 , 
and the fitted / (1710) is only 0.43 ± 0.21% of lutt+tt-. In the uK+K' data 
presented here, the /o(1710) intensity is (38 ± 6)% of the data within the 
same acceptance as for c<J7r + 7r~; the error is almost entirely systematic, and 
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Fig. 3. (a) The projection on to M(K + K ) from an alternative fit using /2(1710), 
(b) the contribution from J p = 2 + . 

covers all alternative parametrisations of the a amplitude and removing the 
Ki(1400). The branching fraction for J/ip -> w/ (1710), / (1710) -> K + K~ 
is (6.6 ± 1.3) x 10~ 4 . We find at the 95% confidence level 



Si2(/o(1710) 



7T7T 



BR(f (17lO) - KK) 



< 0.11, 



where all charge states for decay are taken into account. 

One caveat is necessary. In our study of J ftp — > <J)tt + tt~ and (f)K + K~ [7], defi- 
nite evidence is found for an / (1770), distinct from / (1710) and decaying to 
mi (and possibly weakly to KK). There is a remote possibility that /o(1710) 
and /o(1770) are both present in utttt data but cancel by destructive interfer- 
ence. Such a cancellation would require that they have the same magnitudes 
but opposite phases. Even then, the cancellation is incomplete, because they 
have different masses and widths. Allowing for this possible cancellation, the 
upper limit of the branching ratio given in eqn. (1) could increase to 0.16 if 
the magnitudes happen to be equal, which is unlikely. 

The peak in Fig. 2(c) at ~ 1550 MeV/c 2 may be fitted with either ^(1525) 
or / 2 (1565), or both. Spin 2 is required by non-isotropic decay angular dis- 
tributions; a fit with an /o with the same mass and width gives a worse log 
likelihood by 64. Also no / (1500) is visible in the ojix+tx~ data of Ref. [1]. 
If the peak is fitted with /2(1525), the branching fraction is close to that for 
^(1270) — > KK. However, because of interferences between helicity ampli- 
tudes, the branching fraction could be a factor 2 larger or smaller. If the peak 
is fitted with / 2 (1565), the branching fraction is similar to that of /2(1565) 
in UJ7T7T data, but again could be a factor 2 larger or smaller. The fit shown 
in Fig. 2 uses /2(1525). The branching ratio of ^(1270 ) between KK and txtx 
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is (5.2 ±2.5)%, consistent with the range of values quoted by the Particle 
Data Group [9]; again the error arises from flexibility in interferences between 
helicity amplitudes. 

There is a further feature at ~ 2150 MeV/c 2 in the K + K~ mass spectrum. 
Some spin > 2 component is required by non-isotropic decay angular distri- 
butions. An optimum fit to present data may be achieved with a mass of 
2150 ± 20 MeV/c 2 and a width T = 150 ± 30 MeV/c 2 ; these values are within 
a few MeV/c 2 of PDG values. Errors are mostly statistical but also cover 
changes when the small amplitudes are omitted from the fit. The data do not 
rule out the possibility of spin 4, but the fit is consistent with the known 
/ 2 (2150) [9]. 
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Fig. 4. Angular distributions for Mkk > 2000 MeV/c 2 for angles x> &u>i and f3 w 
defined in the text; histograms show the fit and the lower shaded histograms the 
background, taken from sidebands. The dashed histograms show the acceptance. 

Fig. 4 shows distributions for four angles after selecting Mkk > 2000 MeV/c 2 . 
The angle x is f ne angle between the decay plane of u — > 7r + 7r~7r° and the 
decay plane X — > KK; Q u is the production angle of the u in the J/ip rest 
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frame. The angle olk is the decay angle of the K in the rest frame of X, 
taken with respect to the direction of the recoil u>; (3 W is the angle between 
the normal to the u decay plane and the beam direction. The distribution for 
cos ax is distinctly non-isotropic, although after integrating over all but one 
of the angles, much of the spin information is lost; the full amplitude analysis 
is much more reliable than projections on to individual angles. The dashed 
curves illustrate the acceptance. The shaded histograms at the bottom of each 
panel show background, which is taken from sidebins. 

A marginal improvement of 21 in log likelihood may be obtained by adding 
/o(2100) — > K + K~ . However, this is not sufficient to be sure of its presence, 
so it is omitted. 

~v 400 
^ 300 

<2 
s 

£ 200 
100 



1.5 2 2.5 

M((oK) (GeV/c 2 ) 

Fig. 5. The projection on to M(ljK ) omitting the loK signal at 1945 MeV/c 2 . 

The ujK mass distribution is fitted poorly unless a component decaying to 
tuK is included with M ~ 1945 MeV/c 2 , T ~ 270 MeV/c 2 . Fig. 5 shows the 
poor fit without this additional amplitude; no background appears capable 
of explaining this effect. The optimum fit is obtained with orbital angular 
momentum I = in the ujK system, i.e. J p = 1 + ; this improves log likelihood 
by 113, and contributes 9.4% of all events. The observed isotropic decay can be 
fitted not only by J p = 1 + , but through conspiracy between several production 
amplitudes for J p = 2" and 0~. The known i^2(1820) with J p = 2~ and 
.ff (1830) with J p = 0~ [9] do not alone give an adequate fit but may make 
some contribution. Our conclusion is that some ujK contribution is needed in 
this mass range, but cannot be identified cleanly and could be a superposition 
of more than one resonance with J p = 1 + , 0~ and 2~. Conclusions about fo 
and /2 components are insensitive to this ambiguity. At lower masses, inclusion 
of Zfi(1400) — > uoK also gives a significant improvement of 81 in log likelihood. 

In summary, the main features of the data are peaks which may be attributed 
to / (1710), / 2 (2150), / 2 (1270) and either $(1525) or / 2 (1565). An upper limit 
of 0.11 is set on the ratio BR{f (1710) -> 7r7r]/ 5 ^[/o(1710) -> KK]. This 
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upper limit could rise to 0.16 if there is a fortuitous cancellation of /o(1710) 
and /o(1770) in ujtttt data in both magnitude and phase. 
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